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^ '• Abstract 

The energy spectrum and corresponding wave functions of a flat quantum 
f^ ' dot with elhptic symmetry are obtained exactly. A detailed study is made 

of the effect of ellipticity on the energy levels and the corresponding wave 
CN . functions. The analytical behavior of the energy levels in certain limiting 

^ I , cases is obtained. 
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I. INTRODUCTION 



Quantum dots have been the subject of both experimental and theoretical research in 
^ . recent years [|I|. They have been successfully created experimentally by applying lithographic 

O '. and etching techniques to impose a lateral structure onto an otherwise two-dimensional 

electron system. These structures introduce electrostatic potentials in the plane of the two- 
dimensional electron gas, which confine the electrons to a dot region. The energy levels of 
^ ' electrons in such a quantum dot are fully quantized like in an atom, and therefore they are 

also referred to as artificial atoms. For a theoretical description of such systems it is essential 
to know the confinement potential. Often a parabolic potential is used. Another simple 
model, which also has the advantage of representing a finite confinement, is a piecewise 
constant potential with a circular boundary 0. This type of potential represents more 
closely the situation in self-assembled quantum dots |I[]. The latter model is used, but 
extended to an asymmetric situation, i.e. an elliptic dot. 

There have been many calculations of energy levels in rectangular dots and wires, see 
for example ||^-§[, but in general, the problem of even one particle confined in such a po- 
tential profile is not exactly solvable, because the corresponding Schrodinger equation is not 
separable. Nevertheless one usually assumes separability, although this approximation is 
only acceptable for deep wells, where only a negligible part of the wave function is situated 
outside the dot. In our study we use an elliptically shaped potential profile for which exact 
solutions can be found, even when we take ^| interface 7^ 0, or a potential Vq ^ oo outside 
the dot. We are especially interested in the results for elongated dots, which could then be 



compared with the case of a rectangular slab or even a quantum wire. Our results are also 
applicable to elliptic quantum wires. 

In the literature one can find many references to calculations on elliptic problems, we 
found one particularly useful @], in which the diffraction of waves by ribbons and by slits is 
studied exactly. The problem of a three-dimensional ellipsoidal (3D) quantum well has been 
solved exactly in Refs. [^ and [§ for ^I'linterfacc ~ 0- '^° '^^^ knowledge, the corresponding 
2D problem for finite potential barriers has not been solved. 

The paper is organized as follows. In Sec. II, we obtain the analytical solution for 
the energy and the wave functions. The numerical results are presented in Sec. III. Our 
conclusions are given in Sec. IV. A representation for the Mathieu functions is listed in the 
appendix. 

II. THE PROBLEM AND EQUATIONS 

A particle with mass m is confined in a two-dimensional hard- wall potential well of finite 
height with elliptic shape as shown in Fig. 0. The potential profile we consider is constant 
everywhere, V{x,y) = inside the dot (i.e. x^ /a? + y'^ jh^ < 1) and V{x,y) = Vq outside the 
dot. Stationary energy levels and wave functions are found by solving the time-independent 
Schrodinger equation 

- ^^'^(^' y) + ^(^' ^)^(^' y) = ^^(^' y)- (i) 

Transformation from Cartesian to elliptical coordinates x = / cosh u cos f, < n < oo, 
y = f sinhusmv, < v < 2tc, where / = \faF^^^\? is the focal length, enables us to write 
the boundary of the ellipse in a one- variable equation u = U = arctanh (b/a). Introducing 
the electron momentum k"^ = ^{E — V) the equation to be solved becomes 

irT + irT + ^cosh2M^-^cos2t;^ = 0. 2 

ou^ ov^ 2 2 

Proposing a solution \I'(ti,f) = F{u)G{v) we find that the problem is separable in two 

equations in u and v, coupled by a separation constant c and the variable k. Denoting 

q = pk'^/4: we obtain, 

, o - (c - 2g cosh 2u)F(u) = 0, (3a) 

, \ ' +(c-2qcos2v)G(v) = 0, 3b 

which are the characteristic equations for the Mathieu functions. Physical solutions to these 



equations are well-known [|T4| 



G^{v,q) = \''f^'\ (4a) 

{sei{v,q), ' 

fr(«.«) = ^''1!;!"'*'' (4b) 

Ms;' ^(M,gi), 

Frin,^ = CS^'-^: (4c) 

MS;' '(M, -92), 



where / = 0, 1, 2, . . ..The normahzation factors are discarded and will be incorporated later in 
the total wave function \E'(-u, v). The Mathieu functions cei{v, q) and se;(f , q) are respectively 
even and odd in v and are the analogue of the cosine and sine functions in the angular 
solution of the circular dot problem. For the -u-problem we make a distinction between 
solutions inside and outside the dot. The former show oscillatory behavior and are the 
analogue of the Bessel functions Ji{p) in the radial solution of the circular problem. The 
latter are monotonically decaying functions and are the analogue of the modified Bessel 
functions Ki{p) in the circular problem. In choosing these particular functions we used the 
boundary conditions that the wave functions need to be finite when u = and u ^ oo for 
obvious normalization reasons. Outside the dot E < Vq and the variable q becomes negative, 
which we noted explicitly. In the appendix explicit expressions for the Mathieu functions 
are given. 

The wave function \1/(m, w) and its derivatives must be continuous everywhere. This 
condition requires that at every point on the border (U, v) of the quantum dot we have 



dFriu,q)/dul^^^^^^^ _ dFr{u,q)/dul^^^^_^^ 
FriU,q,) FriU^-q^) 

where the -u-independent functions Gi(y,q) cancel out. 

We introduce a = y/ab, which is related to the area of the ellipse and define 



(5) 



2m 
^ = aki^ 1] = crk2, 1 = (y\l -ryK) , (6) 

with kl = ^E and A;| = ^(Vq — E). The energy levels will be expressed in units of Vq. 

Clearly ^^ + r]'^ = 7^ and gi = ^^j^ = 1^2^^ = |(f — -)C^- We are therefore left with two 
equations for even and odd solutions in one variable ^ to be solved numerically 



9Mc|^^ /du {u, i(f - !)a 5Mcr^ /du {u, -i(f - !)(7^ - e)) 

u=U 



u=U 



Mcr^ {u, i(f - ^je) Mc}^^ {u, -i(f - ^)(72 - e)) 



9Ms|^^ /du {u, i(f - !)a 5Msr^ /du {u, -i(f - !)(7^ - e)) 

u=U 



u=U 



(7a) 



(7b) 



Msp^ {u, Kf - ^je) Msf) {u, -i(f - !)(7^ - e)) 

Solutions are indexed ^n,h n = l,2,...;/ = 0,l,.... The energy spectrum is then given by 

E^i = —kl =^Vo. (8) 

In the limit of an infinitely high barrier, Vq ^ 00, the boundary condition will simply 
become -F/'^('U = f/, gi) = 0, or equivalently, 

la 6^ 



^-P{u,]CT--)e] =0. (9b) 



When the particle has a different effective mass inside (mi) and outside the dot (7712) the 
boundary condition Eq. (^) has to be modified into 

dFr/du{u,q,)l^^ ^ dFr/du{u,-q,)l^^ 

The exphcit expressions for the wave functions are, for the even functions, 

.j,e(„,A_j K,i^c\^\u,qi„,i)cei{v,qi), u<U, 

[ K,i Mc; [u, -g2„,, ) ce; [v, -g2„,, ) , u > U, 

and for the odd functions, 

xi,o(„,A_j ^n,iMsp^(u,gi )sei(w,gi„,), u < U, , . 

[ K.,i MsJ '{u, -q2,J sei{v, -q2,J, u > U, 

where N^ ^ and N° ^ are determined by normalization of the complete wave function. 

III. NUMERICAL RESULTS 

In Fig. H we show the dependence of the energy of the bound states on the strength of 



the well 7 = -1 /^afeVo, which depends both on the potential height Vq and on the area of the 

dot ab. For clarity we have limited all our graphs to / < 7. For a/b = 1 the results coincide 
with our calculations for a circular dot. When the ratio a/b increases, the degeneracy of the 
/ 7^ energy levels is lifted and they split in separate levels. This is a result of a decrease in 
symmetry of the system in the transition from a circular to an elliptic dot. With increasing 
a/b, energy levels corresponding with even states lower, while those corresponding with 
uneven states rise in energy. This behavior is more clearly visible in Fig. ^, where we have 
plotted the energy levels against the ratio a/b, with 7 fixed. For small and narrow dots, only 
the even n = 1 states remain, while the uneven states become unbound. The total number 
of states stays approximately equal. In Fig. ^ it is also clear that even the smallest and 
narrowest dot holds at least one confined state. The energy of this lowest state is shown in 
Fig. H as a function of the eccentricity a/b. There is only a small increase of the energy 
level with increasing eccentricity. 

The asymptotic behavior for the energy of the ground state in a shallow 2D well can be 
obtained analytically. Because of the weak dependence of this energy on the eccentricity 
we assume for simplicity a/b = 1. Outside the well the wavefunction is then given by the 
modified Bessel function, 



^{p) = Ko{J—\Eo\p], (13) 




where \Eo\ = Vq — E is the binding energy in the well and p = (x^ + y^) • For small \Eq\ 
this wave function can to first order be approximated by 




^ip)^-^nr-j^\E,\p], (14) 



with g = 0.57721 . . . the constant of Euler-Mascheroni. Inside the potential well the ground 
state wave function is almost constant. Integrating this function over the Schrodinger equa- 
tion we obtain the binding energy 

with the prefactor 4/e^5 = 1.26096 .... In Figs. ^ and ^ the ground state energy for small 7 
can be approximated by 

E/V, = l-^-^e-''-\ (16) 

The results for an infinitely high confinement potential are shown in Fig. |^. The behavior 
is comparable to the finite well case, except that all the energy levels will eventually rise 
with increasing eccentricity. 

Notice that in Fig. ^ we find that for large eccentricity the energy levels appear in bands 
where the levels are practically equidistant. This can be understood from the following 
simple consideration. If b/a -C 1 we can use an adiabatic approximation and assume that 
the motion along the y-direction is much faster than along x. For |x| < a the electron feels a 
hard wall potential along the y-direction of width W = 26-^/1 — x"^ /a? which has energy levels 



o2 



Ey = -^^^ny. This results into a potential along the x-direction V{x) = f^^TYz^JT^ — 

|^^n^(l + x^/a^) which near its minimum is parabolic. The resulting energy levels become 

En^,ny = '^^^1 + ^^^^y {"^x + 1/2)7 which for fixed Uy results into an equidistant set of 
levels as found in Fig. ^ for a/h ^ 1. The dependence on the eccentricity becomes more 
clear if we introduce a = vab, which was taken fixed in Fig. |^, into the above expression, 
which leads to 

Notice that the linear increase with a/b results from the zero point motion along the y- 
direction. 

To see the effect of different well (mi) and barrier {1712) masses we present in Fig. ^ the 
results for the bound states in the case of three different mass ratios fi = m2/mi = 0.5, 1 
and 2. There is a decrease of the energy levels, i.e. increase of binding, when the effective 
mass of the electron outside the dot is higher than inside the dot. The / = levels become 
bound at the same value for 7, while this is no longer the case for / > 0. 

The density distribution of the wave functions corresponding to the different bound 
states, denoted by {n,lY or {n,l)°, are shown in Fig. |^ for a/b = 5 and 7 = 11. The even 
n = 1 states form the lowest levels and they have all their extrema of the wave function 
located along a row in the direction of the longest axis of the ellipse. In this particular case, 
a second row of extrema appears only at the sixth state (1, 1)°, which is odd. When n = 2, 
we see three rows appear. In all cases the number of extrema in a row is determined by I. 
The general pattern is that the lowest levels are those for which the position of the extrema 
of the wave function best fit the shape of the quantum dot. 



IV. CONCLUSIONS 

By expressing the wave functions in the appropriate coordinate system, viz. elhptical 
coordinates, we were able to find analytical solutions for the Schrodinger equation describing 
an elliptical quantum dot with finite height hard walls. These functions are the Mathieu 
functions. The condition of smoothness of the wave function at the boundary of the dot, 
results in an expression for the energy levels. The main effect of the elongation of the dot 
from a circular to an elliptic shape is a lifting of the degeneracy of the / > 1 levels. We found 
exact results for the energy levels and the wave functions for arbitrary quantum numbers. 
The basic solutions found here may be the basis for more advanced problems with elliptic 
geometry, e.g. many particle quantum dots and quantum dots in a magnetic field. The 
results can readily be expanded to the case of an elliptic wire. The particle can then move 
freely along the z-direction of a wire with an elliptic cross- sect ion. For the energy states we 
find 

C:! rtK y, /-I o\ 

En,l,k^ = 7; \- En,l, (18) 

2m 
where En^i are the above discussed levels and the wave functions are 

^{x,y,z) = e'''^'^n,i{x^y)- (19) 
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APPENDIX 

Solutions to Eq. (|3b| ) are found in the form of a Fourier series. Physically relevant 
solutions obey periodic boundary conditions 

G{v) = G{v + n), (A.la) 

G{v) = G{v + 27t). (A.lb) 

Even periodic solutions are expanded in a Fourier series of cosine functions, while odd 
periodic solutions are expanded in a Fourier series of sine functions. In this way, we find the 
four fundamental solutions, which are called the Mathieu functions |]l4|-p!6[] 



ce2i{v,q) = ^ A|(g) cos2jv, (A.2a) 



i=o 



ce2i+i{v, q) = J2 ^I'i'M cos(2j + 1)^;, (A.2b) 

i=o 

6 



se2i+i{v, q) = J2 Bl'j'+M sin(2j + l)v, (A.2c) 



i=o 



^^2i+2{v, q) = Yl BlY+l^q) sin(2j + 2)v. (A.2d) 

i=o 



By substituting these series into Eq. (pE|), we find a set of recursion relations for the 
expansion coefficients, e.g. for A^^i^q) 

cAl' - qAl' = 0, 
ic-A)Al'-qi2Al' + Af)=0, (A.3) 

[c - {2ir] Al - q{Al_, + A%^,) =0 (j > 2). 

Numerical methods to calculate the Mathieu functions are described in detail in the 
literature PHT^,P^. Solutions to Eq. (^) are written as an expansion in Bessel functions. 
With Ml = Y^e~" and U2 = y/qe"^ we find 

^ oo 



oo 



^1 WJ ,=0 



oo 



^ oo 

Msg^2(«, 9) = ;^2lT27T E(-l)'"''^St2(9) [^,(«i)^,+2(«2) - J,+2(«l) J,(«2)] , (A.4d) 

2 *^"' j=0 



and 

2i(-l, 

j=0 



Mcg^K-g) = -^-j^Y.'^'i^^W^^^K.iu^), (A.5a) 



McgJ.i(n, -q) = -^mr^^ E ^StK?) [Ij{ni)K,+i{u2) - I,^,{u,)K,{u2)] , (A.5b) 

4?^i(«, -g) = \^'y^Al+^\iq) [I,{m)K,^,iu2) + I,^,{u,)K,iu,)] , (A.5c) 



Msgi,(«, -g) = -1^^ ^4X\(g) [/,(«i)ir,+i(«2) + /, 



27 f— iV °° 
Msgi2(n,-g) = _-^^^i?g+2(g) [I^{u,)K,^,{u2) - I,+2{n,)K,{u,)]. (A.5d) 

To obtain the derivatives of the Mathieu functions we used these series expansions of the 
Mathieu functions and applied the analytical formulas for the derivatives of Bessel functions 
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FIGURES 

FIG. 1. The potential profile V{x,y) is zero inside the elliptic dot and equal to Vq outside the 
dot. The long axis of the ellipse is 2a, the short axis is 26. 

FIG. 2. The confined even and odd states (n, iy'° of an elliptic quantum dot with eccentricity 
a/b as a function of the strength of the well 7. All dots have at least one bound state, (a) For 
a/b = 1 the even and odd states are degenerate, (b) When a / 6 this degeneracy is lifted for I > 0. 

FIG. 3. The confined energy states (n, 1^'° of an elliptic quantum dot with strength 7 = 9 as 
a function of the eccentricity a/b. All the even n = 1 states decrease with increasing a/b. 

FIG. 4. The energy of the ground state (1,0)"^ of shallow elliptic quantum dots with small 
strengths 7 = 1, 1.5, 2, 3 and 4 as a function of the eccentricity a/b. 

FIG. 5. The energy states in an infinite elliptic quantum well as a function of the eccentric- 
ity a/b. The effective electron mass was taken to be m* = 0.041 rrie and the area of the dot 
a = {ab) ' = 10 nm. 

FIG. 6. The confined energy states of an elliptic quantum dot with eccentricity a/b = 5 as a 
function of the strength of the well 7 for different mass ratio /i = 1712 /mi where nii (771,2) is the 
effective mass of the electron inside (outside) the dot. For clarity we have split up the graph for 
different /. 

FIG. 7. The probablity density for the eigenfunctions {n,lY'" of all the bound states in an 
elliptic quantum dot with strength 7 = 11 and ratio a/b = 5 in order of increasing energy. The 
lowest states have 77 = 1 and are all even: the maxima in the probability are along the long axis of 
the ellipse. Only higher energy states can have more rows of maxima. 
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